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ABSTRACT 

We present evidence for a Spitzer-selected luminous infrared galaxy (LIRG) behind the Bullet Cluster 
(1E0657-56). The galaxy, originally identified as a multiply imaged source using IRAC photometry, 
has a spectral energy distribution consistent with a highly extincted (Ay ~ 3. 3), strongly star- fo rming 
galaxy at z — 2.7. Using our strong gravitational lensing model presented in Bra dac et alj (|2006[ ). we 
find that the magnifications are « 10 — 50 for the three images of the galaxy. The brightest and 
faintest images differ by a factor of 3.2 in magnification. The implied infrared luminosity is consistent 
with the galaxy being a LIRG, with a stellar mass of M* ~2x 10 10 M© and a star formation rate of 
~ 90 M© yr" 1 . With lensed fluxes at 24/im of 0.58 mJy and 0.39 mJy in the two brightest images, 
this galaxy presents a unique opportunity for detailed study of an obscured starburst with a star 
formation rate comparable to that of L* galaxies at z > 2. 

Subject headings: galaxies: evolution, starburst — gravitational lensing - galaxies: clusters: general 



1. INTRODUCTION 

Measurements of both the star formation history of 
the universe and corollary build-up of stellar mass 
have established t hat the star formation rate peaks at 
1 < z < 3 (eg. iMadau et al.l H9M iLillv et al.l IIMI 
Dickinson et al.ll2003l : iRudnick et alj l2003UReddv et all 
2008L IWilkins et alJ l2008f T Spitzer 24/xm observa- 
tions further indicate that star formation at this epoch 
is dominated by luminous and u ltraluminous infrared 
galax i es (LIRGs and ULI RGs; iPerez-Gonzalez et alJ 
l2005t iLe Floc'h etafl 120051 ). ULIRGs are sufficiently 
bright to facilitate spectroscopy and detailed analyses 
(iDaddi et alj|2005l lYan et al.ll20ol IValiante et all 120071 : 
iPope et al.ll2008h : however, these galaxies rep resent only 
the most massive tail o f the galaxy population (iDev et all 
[2^lD^e~eTaDl2^Ml . 

In contrast, LIRGs have properties more similar to 
the overall galaxy population, with stellar masses and 
star formation rates comparable to th ose seen for UV- 
selected star-formin g gala xies at z ~ 2 (|Reddv fc Steidell 
120041 iReddv et all 12006ft . The intrinsic faintness of 
LIRGS however precludes both high-fidelity multiwave- 

1 This paper includes data gathered with the 6.5 meter Magel- 
lan Telescopes located at Las Campanas Observatory, Chile, the 
Hubble Space Telescope, and the Spitzer Space Telescope. 
Electronic address: anthony@astro.ufl.edu] 



length photometry and spectroscopic programs at optical 
and infrared wavelengths at z > 2. 

Strong gravitational lensing enables observations of in- 
trinsically fainter galaxies than is otherwise possible, and 
several recent programs have begun to exploit lensing 
by galaxy clusters to probe the proper ties of infrared 
and submillimeter lumin ous galaxies (e.g. iKnudsen et al.l 
120081 iRigbv et al.ll2008i and references therein). The 
main limitation of this approach is simply the small num- 
ber of known lensed galaxies that are luminous at these 
wavelengths. 

In this paper we present evidence for a strongly lensed, 
luminous infrared galaxy that is triply imaged by the 
Bullet Cluster. This galaxy is the only strongly lensed 
source for which the initial detection was made with 
Spitzer at mid-infrared wavelengths, and due to its large 
magnification provides a window onto the properties 
of lower luminosity infrared galaxies than have previ- 
ously been studied in this redshift regime. In previ- 
ous papers our team has explored the physical prop- 
erties of the Bullet Cluster, 1E0657-56, measuring its 
matter distribution and the properties of the X-ray gas 
dMarkevitch et alJ^OOl [2001 iClowe et alJ l2(¥)i l200l 
iBradac et al.l 1200*6 ). and const raining the dark mat- 
ter self-interaction c ross-section (Markev itch et aDl2004J : 
iRandall et alJ [20071 ) . The object that is the subject of 
this paper was first identified as a doubly lensed source 
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Table 1. Observed Fluxes and Magnitudes for Lensed Images 
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llux^Jy) 
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B/C 


F606W a 


<0.29 


>25.26 


<0.29 


>25.26 


<0.29 


>25.26 






F775W 


<0.24 


>24.67 


<0.24 


>24.67 


<0.24 


>24.67 






F850LP 


<0.41 


>24.11 


<0.41 


>24.11 


<0.41 


>24.11 






Jc 


<1.30 


>23.62 


<1.30 


>23.62 


<1.30 


>23.62 






K s 


<5.35 


>22.08 


<5.35 


>22.08 


<5.35 


>22.08 






3.6/im b 


13.6 ±0.6 


21.06 ±0.05 


21.0 ± 1.5 


20.59 ± 0.08 


7.4 ± 1.7 


21.73 ±0.25 


1.54 ±0.13 


2.84 ±0.68 


4.5/im 


23.4 ±0.9 


20.47 ±0.04 


32.7 ± 1.5 


20.11 ±0.05 


10.6 ± 1.7 


21.34 ±0.18 


1.40 ±0.08 


3.08 ± 0.52 


5.8/im 


38.4 ± 1.7 


19.94 ±0.05 


57.9 ± 1.9 


19.49 ± 0.03 


16.8 ± 1.7 


20.84 ±0.11 


1.51 ±0.08 


3.45 ± 0.37 




46.4 ±2.5 


19.73 ± 0.06 


67.1 ± 2.5 


19.33 ± 0.04 


20.6 ±2.5 


20.62 ± 0.13 


1.45 ±0.09 


3.25 ± 0.41 


24^m 


390 ± 20 


17.42 ±0.05 


575 ± 20 


17.00 ± .04 


175 ± 20 


18.29 ±0.12 


1.47 ±0.09 


3.29 ±0.39 



in one of these papers (jBradac et al.ll2006l §6 and Fig. 
8), and is independen tly detected as a millimeter source 
(jWilson et al.l l2008af ) . Here we incorporate new HST, 
Spitzer, and Magellan data, in a detailed analysis of this 
object. The data are presented in §2 and are used in $3] 
to estimate the redshift, magnification, stellar mass, and 
star formation rate. In this section we also present evi- 
dence for a newly discovered, third image of this galaxy. 
We summarize our results in §31 

2. PHOTOMETRIC DATA AND MEASUREMENTS 

We use the combination of Sv itzer IRAC (iFazio et al.l 
2004) and MI PS (|Rieke et all I2004D . 
HST ACS (|Ford et al.ll200l . and ground-based near- 
infrared (NIR) observations to constrain the spectral 
energy distribution of the multiply imaged source. The 
subsections below describe the data and photometric 
analysis associated with each facility. 

2.1. Spitzer IRAC 

We originally detected this source as a doubly imaged 
object in Spitzer IRAC data obtained on December 17- 
18, 2004. The data from this program include imaging in 
all four IRAC bands (3.6/im, 4.5/im, 5.8/im, and 8/im). 
These data were taken with a cycling dither pattern with 
medium scale factor and 100s frame time during an 8720s 
duration Astronomical Observation Request (AOR). The 
effective exposure times are 4 ks in each filter. 

We process the data using MOPEX (jMakovoz et all 
2006), with a final pixel scale of 0'/86. Before measuring 
aperture fluxes, we first use GALFIT (jPeng et al.ll2002f ) 
to model and subtract a cluster elliptical that lies directly 
between the two lensed images (object #1 in Figure [l}, 
using a nearby, isolated star as the input PSF for GAL- 
FIT. The structural parameters for the galaxy are held 
fixed to values derived using the ACS data (see below), 
with only the position and magnitude permitted to vary. 4 
There is also a star between the two lensed images (ob- 
ject #2 in Fig. [T|) that is detected at 3.6/im and 4.5/im, 
but is fainter at these wavelengths than the lensed im- 
ages. Due to its faintness, we mask this star rather than 
model it with GALFIT. We also mask two other nearby 
sources that lie within the background apertures (#3 and 
#4 in Fig. [J). 

4 The position is permitted to vary at the subpixel level to mini- 
mize errors in the subtraction due to residual mis-registratioin be- 
tween the images. Fixing the position does not qualitatively alter 
our results. 



We then perform aperture photometry using an aper- 
ture of radius 2.4", with a background annulus extend- 
ing from 2.4 — 7.2", applying the point source aperture 
corrections given in the IRAC Data Handbook. 5 We 
note that the total magnitude may be underestimated 
by ~ 10% due to the modest spatial extension of the 
source in the IRAC data, but apply no additional cor- 
rection for this factor. We measure the flux in an ensem- 
ble of off-source apertures to compute the photometric 
uncertainties. The resultant photometry is given in Ta- 
ble Q] for these two images, which we denote as A and B 
(Fig. [Hand [2]). Image B is brighter than A by a factor 
of 1.5, with consistent flux ratios in all bands (Fig. [3]), 
as required for a multiply imaged source. The coordi- 
nates are (a^ooo, ^2000) =(06:58:38.0,-55:57:02) for image 
A and (02000, <^20oo)=(06:58:37. 1,-55:57:06) for image B. 

In Table [1] we also present photometry for a newly dis- 
covered third image of the galaxy, hereafter denoted as 
image C (see Fig [IJ , which is discussed in greater detail 
in §3.21 The IRAC photometry for image C is obtained in 
the same fashion as for images A and B, in this case mod- 
elling and subtracting two nearby stars and one nearby 
galaxy (objects #5-7 in Fig. 1). 

2.2. Spitzer MIPS 

The MIPS 24/im data were acquired on November 30, 
2007. Observations were taken in small scale photometry 
mode using a 3 x 3 raster map with each position offset by 
half the array. The frame time for individual exposures 
is 30s. 

We process the data using MOPEX, with a final 
pixel scale of 1"22. At this wavelength images A, B, 
and C are all clearly detected. We use APEX (Astro- 
nomical Point source Ex tractor; iMakovoz et al.l 120021 : 
iMakovoz &: Marleaull2005f ) to perform point source pho- 
tometry, deblending images A and B. The fluxes of all 
three lensed images, which are derived via PSF fitting, 
are reported in Table [TJ The quoted uncertainties in- 
clude uncertainty associated with background subtrac- 
tion, which dominate over the statistical uncertainties 
reported by APEX. The MIPS point-source photome- 
try does not suffer from contamination by foreground 
sources. The spectral energy distributions of cluster el- 
lipticals fall rapidly between 8/im and 24/im and the 
nearby sources modelled at shorter wavelengths are not 

5 See http://ssc.spitzer.caltech.edu/irac/dh/ A 2.4" 
aperture is the smallest tor which corrections are given. 



A Multiply Imaged LIRG Behind the Bullet Cluster 



3 




Fig. 1. — Top- An image showing the location of image C relative to images A and B. In this image the galaxy between images A 
and B has been subtracted for clarity. The field of view is 65" X 50". Bottom - F850LP cutouts of the regions surrounding lensed images 
A and B (left) and image C (right). The crosses denote the locations of each image; the solid curve in the left panel is the critical curve 
from the z = 2.7 magnification map. The objects detected in the F850LP image that lie closest to the IRAC coordinates for images A and 
B are numbered 1-4. Object 1 is a cluster elliptical galaxy and object 2 is a star. The two fainter galaxies (3 and 4) are offset from the 
IRAC detections by 0.8" and 1.5", whereas the relative astrometry is good to 0.25", and can thus be excluded as optical counterparts to 
the lenses. The field of view is 17" X 12" in both panels. For all images north is up and east is to the left. 



detected at 24/im. 

2.3. HST ACS 

We observed the Bullet Cluster with the Advanced 
Camera for Surveys (ACS) through the F606W, F775W, 
and F850LP filters. The F606W data were taken on Oc- 
tober 21, 2004, with a total exposure time of 4.7ks. The 
F775W and F850LP data were taken on October 12-13, 
2006, with total exposure times of 10. Iks and 12.7ks, 
respectively. All observations were reprocessed and driz- 
zled to a common coordinate system using custom soft- 



ware (Haggles; Marshall et al. 2008, in prep) based upon 
MultiDrizzle (|Koekemoer et al.ll2002| ). 

For the HST data we perform aperture photometry 
within 1.5" aperture radii. This aperture size is selected 
as a tradeoff between two competing factors. Specifically, 
while smaller apertures yield more stringent lower limits 
on the magnitude, the aperture size must be sufficiently 
large to encompass the total flux from each lensed image. 
Given that IRAC provides the highest resolution data in 
which the source is detected, our information about the 
true physical extent of the lensed images is limited. The 
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Fig. 2. — Cutout images of the region around the lensed galaxy, in order of increasing wavelength starting with the ACS F850LP data. 
In the first two rows we show the images used for the photometry after subtraction of contaminant sources using GALFIT. In the last row 
we show images prior to GALFIT subtraction for a representative sample of wavelengths (F850LP, J,4.5fim, and 8/im). In all panels the 
field of view is 34" X 26". North is up and east is to the left. 



galaxy shows only limited spatial extent at 3.6/zm, indi- 
cating that a 1.5" aperture is sufficiently large to enclose 
the total flux in the ACS imaging. Any future higher res- 
olution detection of this object will enable use of smaller 
apertures and yield improved magnitude limits. As with 
the IRAC data, we determine the photometric uncer- 
tainty using an ensemble of background apertures. 

Before measuring aperture fluxes, we again use GAL- 
FIT to model and subtract contaminant sources near 
images A and B. In the HST imaging this includes ob- 
jects #1 and #2. Near image C there are no objects 
that require subtraction for the HST data. We again 
use a nearby, isolated star as the input PSF for GAL- 
FIT, which has the advantage over Tin yTim 6 of correctly 
repro ducing the red halo in F850LP (jGilliland fc Riesa 
2002), and we recover consistent structural parameters 
for the cluster galaxy in all filters (effective radius r e = 
0.9" and Sersic index n = 4.6). We also mask out galaxy 
=#=3 in the aperture of image B (Fig. [J). This galaxy 
cannot be the optical counterpart to the IRAC detection 
since it offset from the IRAC detection by 0.8", whereas 
the relative astrometry is good to 0.25". Moreover, both 
the consistent flux ratios in all IRAC bands for images A 
and B (see Table Q] and £j3j) and the location of the criti- 
cal curve in the lensing model support the interpretation 
that these are multiple images of the same source. In this 
case, the relative flux ratio of the two images should also 
be preserved in the HST data, and we would detect the 
counterimage in the other aperture at high confidence 
if it were the optical counterpart. For completeness, in 
Table [2] we provide the HST photometry for the objects 
that are labelled in Figure [1] computed using Source Ex- 

6 http: //www. stsci . edu/sof tware/tinytim/t inytim.html 



tractor. 

2.4. Magellan PANIC 

We imaged the centra l region of the Bullet Cluster with 
the PANIC instrument (jMartini et al.ll2004h on Magellan 
on March 06, 2006. Data were obtained in the J c and 
K s filters and photo metrically calibrated t o the 2MASS 
point source catalog (|Skrutskie et al.|[2006f ). with seeing 
of 0.55 — 0.6" in both bands. Similar to the approach 
taken with the other data sets, we use GALFIT to fit and 
subtract off the bright galaxy and star that lie between 
the locations of images A and B in the IRAC data. We 
then measure the flux within the same 1.5" apertures 
employed for the ACS analysis, recovering only upper 
limits at the positions of all three images. 

3. ANALYSIS 

3.1. Spectral Energy Distribution and Photometric 
Redshift 

In Figure [4] we plot the spectral energy distribution 
(SED) for each image of the lensed galaxy. Qualitatively 
the combination of strong upper limits at optical and 
near-infrared wavelengths coupled with IRAC detections 
and a strong MIPS detection argue for the galaxy be- 
ing a dusty starburst at z ~ 2, with the MIPS 24/im 
emission being due to the redshifted PAH features. The 
24/im emission is difficult to explain if z > 3, while the 
galaxy should be detected at NIR or optical wavelengths 
if either the internal extinction is low or the redshift is 
much below 2. 

For a more quantitative answer, we use the pho tomet- 
ric redshift code HyperZ (jBolzonella et all 1200(1 . The 
input spectral templates are obtained using the Char- 
lot & Bruzual 2007 models (jBruzuall 12007) with the 
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Table 2. HST Photometry for Objects Near Lensed Images 



r lg. i 


r OUO Vv 


r / 1 VV 


F850LP 




ID 


(AB) 


(AB) 


(AB) 




Type 


1 


19.89 ± .02 


19.05 ± .03 


18.69 ± 


03 


Galaxy 


2 


19.79 ± .01 


18.95 ± .01 


19.17 ± 


01 


Star 


3 


24.27 ± .07 


23.43 ± .07 


23.39 ± 


06 


Galaxy 


4 


24.42 ± .07 


23.58 ± .08 


24.06 ± 


08 


Galaxy 


5 


20.34 ± .01 


19.50 ± .01 


18.75 ± 


01 


Star 


6 


19.81 ± .01 


18.97 ± .01 


18.60 ± 


01 


Star 


7 


21.27 ± .03 


20.43 ± .03 


20.18 ± 


03 


Galaxy 



Note. — In this Tabic wc quote Source Extractor AUTO magni- 
tudes. The uncertainties arc calculated using artificial stars and galax- 
ies. 



Padov a 1994 evolution ary tracks (jBertelli et all 119941) 
and a iChabrierl ((20031) mass function. The templates 
are denned to have star formation histories identical to 
the d efault synthetic temp lates provided with HyperZ, 7 
and a lCalzetti et alJ (|2000l ) extinction law is employed. 

Since these stellar templates do not include PAH emis- 
sion, for the photometric redshifts we fit only the data 
shortward of 10/im. HyperZ yields a best fit redshift 
z = 2.72+0 ^ for image A (90% confidence; xl = 0.5). 
A small secondary peak in the redshift probability dis- 
tribution is observed at z — 5 (Fig. [5|); however, this 
redshift is implausible because of the high observed flux 
at 24/im. An analysis of image B yields a similar redshift 
z = 2.62l°;22 (xl = 0.97) with no secondary peak. In 
IBradac et alJ (|2006h we speculated that the source might 
be at z > 6, but this possiblity is now excluded at high 
confidence (A^ 2 > 8). The best fit spectral template 
(Fig. [4]) corresponds to a dusty starburst galaxy with 
Ay = 3.3to 8 (90%) and an age of < 30 Myr. As can 
be seen in the figure, the robustness of the photometric 
redshift is largely due to the fact that the IRAC data 
span the 1.6/xm bump at this redshift. 

At the best-fit redshift the 7.7/um and 8.6/im PAH fea- 
tures (blended in this spectrum) are redshifted beyond 
the 24/im window, in which case the observed 24/im emis- 
sion is dominated by the 6. 24/im PAH feature. In a con- 
current program ([Wilson et al.l 12008 af ) have also identi- 
fied this galaxy as a bright millimeter source, and obtain 
a consistent redshift (z — 2.7) via an empirical relation 
for SMGs between redshift and IRAC colors. 

3.2. Magnification and Additional Images 

The magnification map was obtained from the strong 
(information from multiply imaged systems) and weak 
(measuring distortion of background galaxies) gravita- 
tional lensing data. It is the same reconstruction as in 
IBradac et alJ (|2006h . It is performed on a pixelized grid 
and does not assume a specific form of the underlying 
gravitational potential. At the location of image A we 
compute a magnification ~ 25 and at image B we ob- 
tain |/i ~ 50 for a galaxy at z — 1.1 . The critical curve 

7 These templates have exponentially declining star formation 
rates with r = 1,2,3,5,15,30 Gyr for E, SO, Sa, Sb, Sc, and Sd 
galaxies, respectively. There is also a starburst template that cor- 
responds to a single, instantaneou s burst model. For all t emplates 
the metallicity is solar; however, Bolzonclla ct al. (2000) demon- 
strated that the redshift determination is not strongly dependent 
upon metallicity. 
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Fig. 3. — Flux ratios for the three images as a function of wave- 
length. The solid lines correspond to the weighted mean flux ra- 
tios, including data at all wavelengths. The observed flux ratios 
at different wavelengths are consistent to within the photometric 
uncertainties For B/A (triangles) and B/C (circles) these mean 
values are 1.47 ± 0.06 and 3.18 ± 0.23 respectively. The flux ratio 
for A/C, which is not plotted, is 2.16 ± 0.18. 



(points of infinite magnification) passes between images 
A and B (see FigfT]) - their parity is reversed - sup- 
porting the hypothesis that A and B are indeed multiple 
images of the same source. These magnifications change 
by <20% within the range of allowed photometric red- 
shifts (90% confidence interval, see N3.ip . The measured 
flux ratio of the two images is 1.47 ±0.05 which is simlar 
to the flux ratio of ~ 2 given by our lens model. In the 
subsequent discussion we will quote the stellar mass and 
star formation rate in terms of (/ia/25) -1 to reflect the 
inherent magnification uncertainty. 

Given both lensing model and positions of images 
A and B, we can search for additional images of 
this source. We identified image C (Figure [1} at 
(«2000j ^2000 ) =(06:58:33.4,-55:57:29) using the same ini- 
tial lens model, which does not include information from 
images A and B of this object. The photometry for im- 
age C, described in indicates that this image is a 
factor of 2.2 fainter than image A, while our lens model 
predicts a factor of four. Given the uncertainties these 
images are still consistent with being multiple images. 
The best- fit photometric redshift derived for this source 
is z = 2.82tjj;& 

3.3. Stellar Mass 

To esti mate the stellar mass we use the code 
kcorrect (Bla nton fc Roweisl 12007ft. which fits a lin- 
ear combination of Bruzual & Chariot (2003) templates 
based upon Padova 1994 isochrones and spanning a range 
in metallicity (0.005 — 2.5 times solar) and age (1 Myr 
to 13.75 Gyr). The required inputs are the photometric 
redshift from HyperZ and the IRAC photometry, red- 
dening corrected us ing the Ay from HyperZ and the 
ICalzetti et alJ (|2000T ) reddening law. In this analysis and 
^3.51 we focus upon image A, but note that equivalent 
results hold for image B. For z — 2.72 and Ay = 3.3 we 
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1 10 10 2 10 3 

Fig. 4. — Spectral energy distribution for the galaxy. The solid points correspond to the observed fluxes for images A (circles), B 
(squares), and C (triangles), while the open arrows correspond to the 5cr upper limits at optical and NIR wavelengths. The solid circles 
at 1.1 mm are the AzTEC data from Wilson et al. (2008), where we have split their observed flux using the flux ratio derived for images 
A and B (see Fig. [3j) . Horizontal error bars on the 24/jm d ata points denote the width of the filter, and for the AzTEC data correspond 
approximately to the system bandpass (Wilson ct al. 2008b). The solid curve at A < 10/jm is the best-fit s pectrum returned by H yperZ for 
image B, which corresponds to a starburst galaxy. The solid curve redward of lOfim is the template from Chary & Elbaz (2001) that best 
fits the 24/^m flux for image A assuming fiA = 25. This template is redshifted to z = 2.7 and rescaled to image B using the observe flux 
ratio of the two images (Fig. [3J. Note that the AzTEC data, while not included in the fit, is fully consistent with this model. The inset 
zooms in on the wavelength regime covered by Spitzer, showing only the photometry for image B for clarity. 



obtain a stellar mass M* = 1.5 x 10 10 (/m/25) _1 M Q . 8 
Thus, we find that this lensed galax y is massive - simila r 
m mass to LIRGs at lower redshift (jCaputi et al.ll2006f ). 

3.4. Presence of an AGN 

In the sections above we have determined the redshift 
and stellar mass assuming that the observed SED is dom- 
inated by stars. It is true however that dusty starbursts 
and active galactic nuclei ar e difficult to discrim inate at 
the source redshift (e.g., see lBarmbv et al.ll2006l ). 

The simplest discriminator between the two contribu- 
tors is spatial extent - any spatially extended emission 
must be stellar rather than due to an AGN. A visual 
inspection of Figure demonstrates that the IRAC im- 
ages do exhibit a modest extension perpendicular to the 
critical curve, arguing that the flux is not purely from 
an AGN. Next, we consider the Chandra observations to 

8 As noted by Maraston et al. (2006), templates that include the 
contribution of TP-AGB stars to the spectrum can yield stellar 
masses roughly a factor of two lower. These templates are not 
currently implemented in kcorrect. 



search for evidence of AGN activity. We find that the 
lensed source is a non-detection in our 500 ks exposure 
(/ < 3.6 x 10 -16 erg s _1 cm" 2 unabsorbed, 3<r, for 0.5-2 
keV). Comparing with the local X-ray to mid-IR rela- 
tions of lKrabbe et a l. (2001), we find that a local Seyfert 
galaxy of comparable mid-IR luminosity should be more 
than a factor of 10 brighter than this limit, whereas non- 
detection is consistent with the expected relation for a 
starburst galaxy. While neither of the above arguments 
exclude an additional contribution from a central AGN, 
they do argue that we are not looking at a purely AGN 
spectrum. 

Given the 24/im data, we can also consider 
whether this source, which is also a millimeter galaxy 
(jWilson et al.ll2008af ). has mid-IR colors consistent with 
a starburst or AGN. Comparing with the distribu- 
tion of 24/im-8/m i vs. 8/im-4.5/^m colors for SMGs in 
iPope et alJ ((2008), we find that the designation is am- 
biguous. The source li e s clos e to, but outside the regime 
defined in IPope et all |2008) for starburst galaxies, ar- 
guing that the observed SED may be a composite with 
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Fig. 5. — The 68% confidence intervals from HyperZ in the 
redshift-extinction and redshift-age planes. The solid curves denote 
the confidence intervals for image A; the dashed curves are for 
image B. The crosses (plus signs) denote the values corresponding 
to the minimum \ 2 f° r image A (B). There is a small secondary 
peak in the solutions for image A at z = 5. 



AGN and starburst contributions. 

3.5. Star Formation Rate 

The only means of estimating the star formation rate 
with the existing data is via the strength of the PAH 
emission. There are two main caveats to this approach. 
First, any AGN contribution will bias our estimate of the 
star formation rate. Second, there exists large scatter in 
the relation between 8/xm emission and s tar formation as 
traced by other methods (jCalzettil 12001 ). 

Keeping the above caveats in mind, we cautiously pro- 
ceed with deriving a rough estimate of the star formation 
rate. To do so, we first convert the o bserved 24 / zm lu - 
minosity to Lm and then use the local iKennicuttl (1998) 
relation to convert Lir to st ar formation ra t e. Us ing 
the templates and code from iCharv fc"lj ]lbaz ( 200]]) to 
fit the 24/im flux, assuming z = 2.7, we derive a best-fit 
Lib. = 5x 10 11 (^ J 4/25)~ 1 L Q . The corresponding implied 
star formation rate is SFR ~ 9O(^a/25) _1 M yr" 1 . As 
a consistency check, we also derive the rest-frame 8/im 
luminosity and convert to Lm usi ng -LsumZ-kiR = 10 , 

(120081 ) 



consistent with recent results from iRigbv et al 



based upon a combination of local and z ~ 2 galaxies. 
To derive the rest-frame 8/im we use the same spec tral in- 
dex correction (a = 2.296) as in lDev et ail (|2008D . This 
approach yields a qualitatively consistent total luminos- 
ity, %~3x lO 11 ^^)- 1 L Q . 

In F igure H] we overlay the best-fit ICharv fc Elbaz 
(2001) template for image B, redshifted to z = 2.7, at 
wavelengths redward of 10/im. If the magnification is a 
factor of two lower than our canonical value, which boot- 
strap simulations indicate is the maximum by which we 



may be overestimating fiA, this galaxy would lie at the 
borderline betwen LIRG and ULIRG luminosity. From 
their AzTEC millimeter data, IWilson et ail (|2008al ) also 
estimate Ljb, — 10 11 — 10 12 L© for this source. The 
estimated specific star formation rate for this galaxy, 
SFR w 5 Gyr -1 , is comparable t o that of similar m ass 
BM/BX galaxies at z = 1.5 - 2.6 (jReddv et al.ll2006f ). 

4. DISCUSSION 

We have presented confirmation observations for a mul- 
tiply imaged source behind the Bullet Cluster, and iden- 
tified a third, previously unknown image. From our mul- 
tiwavelength imaging we argue that the source is most 
consistent with being a dusty (Ay ~ 3.3), strongly star- 
forming galaxy at z ~ 2.7. At this redshift our mass 
model for the cluster core indicates that the galaxy is 
highly magnified (fiA ~ 25), implying a large intrinsic 
stellar mass of M* ~ 2 x 10 10 M©. We estimate a star 
formation rate of ~ 90 M© yr _1 based upon the observed 
flux in the 24/i band, which we assume to be dominated 
by emission from redshifted PAH features. 

The estimated intrinsic infrared luminosity of this 
galaxy (5 x 10 11 L©) qualifies it as a luminous infrared 
galaxy (LIRG), fainter than the ULIRGs typically stud- 
ied at this epoch. The galaxy is also known to be 
an exceptionally b right SMG (13.5 mJy at 1.1 mm, 
IWilson et ail2008a| ). and thus may be an ideal system for 
studying the connection between different classes of in- 
frared sources at lower intrinsic luminosity than has pre- 
viously been possible for this epoch. We anticipate that 
the galaxy will be detected in scheduled deep HST NIC- 
MOS observations, providing information on the spatial 
extent of the lensed images. 

The Bullet Cluster, due to its large lensing cross- 
section, provides an optimal environment in which to 
identify lensed galaxies such as the one presented here. 
However, even in this rare, massive cluster merger, we de- 
tected only one lensed LIRG. This result highlights that 
the prospects are not good for finding large samples of 
z ~ 2, gravitationally-lensed LIRGS. It is for this reason 
that each case must be highlighted and exploited. 
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